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Abstract 



We introduce in this paper a new four-parameter generalized version of the linear failure 
rate (LFR) distribution which is called Beta-linear failure rate (BLFR) distribution. The 
new distribution is quite flexible and can be used effectively in modeling survival data and 
reliability problems. It can have a constant, decreasing, increasing, upside-down bathtub 
(unimodal) and bathtub-shaped failure rate function depending on its parameters. It includes 
] some well-known lifetime distributions as special sub-models. We provide a comprehensive 

account of the mathematical properties of the new distributions. In particular, A closed-form 
expressions for the density, cumulative distribution and hazard rate function of the BLFR 
is given. Also, the rth order moment of this distribution is derived. We discuss maximum 
, likelihood estimation of the unknown parameters of the new model for complete sample and 

l/^ , obtain an expression for Fishers information matrix. In the end, to show the flexibility of 

(N . this distribution and illustrative purposes, an application using a real data set is presented. 

(N 

MSG: 60E05; 62F10; 62P99. 

Keywords: Beta distribution; Hazard function; Linear failure rate distribution; Maximum 
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1 Introduction 

The linear failure rate distribution with parameters a > and 6 > 0, {a + b > 0) which is 
denoted by LFRD(a, 6), has the cumulative distribution function (CDF) 

G (x) = 1 — exp ( —ax — -x^j , x > 0, (1-1) 
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and probability density function 



(x) = {a + bx) exp ( —ax — ^x^ ) , x > 0. (1-2) 
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Note that if 6 = and a ^ 0, then the LFR distribution is reduced to exponential distribution 
with parameter a (Exp(a)), and if a = and b ^ then we can obtain the Rayleigh distribution 
with parameter b (Rayleigh (6)). A basic structural properties of LFRD(a, 6) is that it is the 
distribution of minimum of two independent random variables Xi and X2 having Exp(a) and 
Rayleigh(6) distributions, respectively (Sen and Bhattachrayya, 1995). 

If G denotes the CDF of a random variable then a generalized class of distributions can be 
defined by 

F{x) = lGi,)ia,f5) = ^^^J^ \''-\l-tf-Ut, (1.3) 

for a > and P > 0, where Iy{a,P) = ^^1^^'^^ is the incomplete beta function ratio and 
By {a, 13) = Jq t'^~^{l — t)^~^ is the incomplete beta function. 

Many authors considered various forms of G and studied their properties: Eugene et al. 
(2002) (Beta Normal distribution), Nadarajah and Kotz (2004) (Beta Gumbel distribution), 
Nadarajah and Gupta (2004) and Barreto-Souza ct al. (2011) (Beta Frechet distribution), 
Famoye et al. (2005), Lcc ct al. (2007) and Cordeiro et al. (2008) (Beta Wcibull distribution), 
Nadarajah and Kotz (2006) (Beta Exponential distribution), Akinsete et al. (2008) (Beta Pareto 
distribution), Silva et al. (2010) (Beta Modified Weibull distribution), Mahmoudi (2011) (Beta 
generalized Pareto distribution), Cordeiro et al. (2011) (Beta-exponentiated Weibull distribu- 
tion), Cordeiro et al. (2011) (Beta- Weibull geometric distribution), Singla et al. (2012) (Beta 
generalized Weibull distribution), Cordeiro et al. (2012) (Beta generalized gamma distribution) 
and Cordeiro et al. (2012) (Beta generalized normal distribution). 

In this article, we propose a new four parameters distribution, referred to as the BLFR 
distribution, which contains as special sub- models: the Beta exponential (BE), Beta Rayleigh 
(BR), generalized linear failure rate (GLFR) and linear failure rate (LFR) distributions, among 
others. The main reasons for introducing BLFR distribution are: (i) The additional parameters 
introduced by the beta generalization is sought as a means to furnish a more flexible distribution, 
(ii) Some modeling phenomenon with non-monotone failure rates such as the bathtub-shaped and 
unimodal failure rates, which are common in reliability and biological studies, take a reasonable 
parametric fit with this distribution, (iii) The BLFR distribution is expected to have immediate 
application in reliability and survival studies, (iv) BLFR distribution shows better fitting, more 
flexible in shape and easier to perform and formula for modeling lifetime data. 

The reminder of the paper is organized as follows: In Section 2, we define the BLFR distri- 
bution and outline some special cases of the distribution. We investigate some properties of the 
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distribution in this Section. Some of these properties are the hmit behavior and shapes of the 
pdf and hazard rate function of the BLFR distribution. Section 3 provides a general expansion 
for the moments of the BLFR distribution. In Section 4, we discuss maximum hkehhood estima- 
tion and calculate the elements of the observed information matrix. Application of the BLFR 
distribution is given in the Section 5. A simulation study is performed in Section 6. Finally, 
Section 7 concludes the paper. 

2 Definition of the BLFR distribution and some special cases 

Consider that g(x) = dG{x)/dx is the density of the baseline distribution. Then the probability 
density function corresponding to (11. 3p can be written in the form 

f = 7^7^^ - ^ (2.1) 

We now introduce the BLFR distribution by taking G{x) in (jl.3p to be the CDF (jl.ip of the 
LFR distribution. Hence, the BLFR density function can be written as 

and we use the notation X ~ BLFR (a, b, a, (3). 

The hazard rate function of BLFR distribution is given by 

h (x) = -— — ^ — — f I _ f —ax — -x^] I exp ( —aBx — — | . (2.3) 

Plots of pdf and hazard rate function of the BLFR distribution for different values of it's 
parameters are given in Fig. [T]and Fig. [21 respectively. 

2.1 Special cases of the BLFR distribution 

1. If /3 = 1, then we get the generalized linear failure rate distribution (GLFR(a, b, a)) which 
is introduced by Sarhan and Kundu (2009). 

2. If /3 = 1 and 6 = 0, then we get the generalized exponential distribution (GE) (Gupta and 
Kundu, 1999). 

3. If /3 = 1 and a = 0, then we get two-parameter Burr X distribution which is introduced 
by Surles and Padgett (2005) and also is known as generalized Rayleigh distribution (GR) 
(Kundu and Raqab, 2005) . 
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4. If a = /3 = 1, then (2.2) reduces to the hnear failure rate distribution (LFR(a,6)) distri- 
bution. 

5. If 6 = 0, then we get the beta exponential distribution (BE(a, q,/3)) which is introduced 
by Nadarajah and Kotz (2006). 

6. If a = 0, then we get the beta Rayleigh distribution (BR(6, /?)) which is defined by Akinsete 
and Lowe (2009) and is a special case of beta Weibull distribution (Famoye et al., 2005). 

7. If the random variable X has BLFR distribution, then the random variable 

y = 1 - exp i-aX - -X^ 



2 

satisfies the beta distribution with parameters a and /3. Therefore, 

b n 

T = aX + -X^ 
2 

satisfies the beta exponential distribution with parameters 1, a and f3 (BE(1, a, /?)). 

8. If a = i and (3 = n — i, where i and n are positive integer values, then the f{x) is the 
density function of ith order statistic of LFR distribution. 

The following result helps in simulating data from the BLFR distribution: If Y follows Beta 
distribution with parameters a and /3, then 



X = G-^ (Y) = < 



-a+y/a^-2blog(l-Y) ^ ^ ^ 



b 

-l^^silzXL if a>0, 6 = 0, 

follows BLFR distribution with parameters a, b, a, and /3. 



For checking the consistency of the simulating data set form BLFR distribution, the his- 
togram for a generated data set with size 100 and the exact BLFR density with parameters 
a = 0.2, b = 0.1, a = 2, and /? = 0.3, are displayed in Fig [3] (left). Also, the empirical 
distribution function and the exact distribution function is given in Fig [3] (right). 

2.2 Properties of the BLFR distribution 

In this section, limiting behavior of pff and hazard rate function of the BLFR distribution and 
their shapes are studied. 

Theorem 1. Let f{x) be the pdf of the BLFR distribution. The limiting behaviour of / for 
different values of its parameters is given bellow: 



i. If a = 1 then lima;_^o f{x) = 

ii. If a > 1 then limjj-^o f{x) = 0. 

iii. If < a < 1 then hma;_^o f{x) = oo. 

iv. Iima;^oo f{x) = 0. 

Proof: The proof of parts (i)-(iii) are obvious. For part (iv), we have 

b ^^2\ ^ ^ n ^ r \ ^ a + hx ( b/3 2 



< - exp \^-ax - 2^ j j <1^0< f{x) < ^ cxp ( -a/3x - —x 
It can be easily shown that 

Urn (a + 6x)exp ( —a^x — ^x"^] = 0. 
X-&00 \ 2 J 

and the proof is completed. ■ 

Theorem 2. Let f{x) be the density function of the BLFR distribution. The mode of / is 
given in the following cases: 

i. If a = 1 and —a + > then f(x) has a unique mode in x = \{—a + ) 

i. If a = 1 and —a + < then /(x) has a unique mode in x = 0. 

ii. If a > 1 then f{x) has at least one mode. 

Proof: The proof is obvious and is omitted. ■ 

Theorem 3. Let h{x) be the hazard rate function of the BLFR distribution. Consider the 
following cases: 

i. If a = 1 and 6 > then BLFR distribution has an increasing hazard rate function. 

ii. If a > 1 and b> then the hazard rate function of the BLFR distribution is an increasing. 

iii. If 6 = BLFR distribution has a decreasing hazard rate function for a < 1(> 1), and h{x) 

is constant for a = 1. 

iv. If a < 1 and 6 > then h{x) is a bathtub-shaped. 



Proof: 

i. If a = 1 then B (a, P) - Bq^x) («, 13) = ^ exp (^-a^x - ^x^j . Therefore 

n [x) = -p X— exp I —apx — —x = {a + ox) , 



^ exp {^—af3x — ^x'^^ 
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which is an increasing and hnear function with respect to x. 
ii. Consider 

bo b/ a\2 



z = ax+-x =^[x + -) -- 
It imphes that z > for x > and also, it is increasing with respect to x. We have x 



^\/lbz + o? — |. Now, rewriting the BLFR density as function of ^ {£) say, we obtain 



Therefore, we have 



and we conclude that the hazard function of BLFR distribution is increasing, 
in. If 6 = then 

log(/(x)) = log (a) — log (i? (a, ^)) + (a — 1) log (1 — exp (— ax) ) — af3x, 

and 

92 (a - 1) exp (-ax) 

&^'°^<^<"» " (l-exp(-»x) f 

Thus we have ^log (/ (x)) > 0(< 0) where a < 1(> 1), which implies the decreasing (increas- 
ing) hazard rate functions in this cases. 

iv. It is difficult to determine analytically the regions corresponding to the upside-down bathtub 
shaped (unimodal) and bathtub-shaped hazard rate functions for the BLFR distribution. How- 
ever, by some graphical analysis wc can shows: bathtub-shaped hazard rate function correspond 
to a < 1 and 6 > 0. the proof is completed. ■ 

3 Some extensions and Moments of the BLFR distribution 

Here, we present some representations of CDF, PDF, and the survival function of BLFR distri- 
bution. The mathematical relation given below will be useful in this section. If ^9 is a positive 
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real non-integer and \z\ < 1, then 



oo 



and if /3 is a positive real integer, then the upper of the this summation stops at /? — 1, where 

(-iyr(/3) 
' r(/3-i)r(j + i)- 

1. We can express (jl.3p as a mixture of distribution function of generalized LFR distributions 
as follows: 

oo oo 

F{x) = J2PjiG{x)r+^ = ^PjG,{x), 

j=0 j=0 

where pj = r r i ^ yp(^j^i'^(^a^j ) and Gj (x) = (G(x))"'^-' is distribution function of a random 
variable which has a generalized LFR distribution with parameters a, b, and a + j. 

2. We can express (|2.2p as a mixture of density function of generalized LFR distributions as 
follows: 

oo oo 

fix) = Y,pM+J)s{x) {G{x)r+^-' = Y^p.g, (x), 

j=0 j=0 

where gj (x) is density function of a random variable which has a generalized LFR distribution 
with parameters a, b, and a + j. 

3. The cdf can be expressed in terms of the hypergeometric function and the incomplete beta 
function ratio (see, Cordeiro and Nadarajah, 2011) in the following way: 

^ = -irr-^ 2Fi{a,l-f3-a + l;G (x)) , 
aB (q, p) 

where ^Fi (a, b; c; z) = EZo / ((c)fcA;!) zK 

4. The kth moment of BLFR distribution can be expressed as a mixture of the kth moment of 
generalized LFR distributions as follows: 

POO POO °° 

E{X'') = / x^f{x)dx= / x'^ Vpi(a+j)g(a;)(G(x))"+^-^(ix 
io Jo ^-^0 

oo oo 

= Vp,- / x^g^{x)dx = y^p,E{X^), (3.1) 

where gj (x) is density function of a random variable Xj which has a generalized LFR distribution 
with parameters a, 6, and a + j. 
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4 Estimation and inference 



Consider Xi, . . . Xn is a random sample from BLFR distribution. The log-likelihood function 
for the vector of parameters = (a, 6, a, /?) can be written as 



£(6>) = ^log(a + 62;0 -nlog(r(a)) -nlog(r(/3)) 

i=l 

n n 

+nlog {T{a + P)) + {a-l)Y, log (1 " exp (U) ) + ^ J] t^, 



(4.1) 

i=l i=l 

where tj = f. The log- likelihood can be maximized either directly or by solving the 

nonlinear likelihood equations obtained by differentiating (14. ip . The components of the score 
vector U (6) are given by 

i=l 



exp {ti 



d 



Xi 



1=1 

d 



i=l ^ ^ '■' i=l 

2, 



(g - 1) x|exp {U) (3 2 
2 l.l-exp(t.) -2^.-- 



C/a (f) = (0) = -ni; (q) + nV' (a + /3) + ^ log (1 - exp {U)) , 



i=l 
n 



Up (6) = —i (0) = -n^ {(3) + nV' (a + /3) + ^ U. 



1=1 



where ip (.) is the digamma function. 

For interval estimation and hypothesis tests on the model parameters, we require the observed 
information matrix. The 4x4 unit observed information matrix J = J{9) is obtained as 



J 



'Jaa 


Jab 


Jaa 


Jal3 


Jba 


Jbb 


Jba 


Jbl3 


Jaa 


Jab 


Jaa 


Jal3 


Jpa 


Jpb 


Jpa 


Jf}fi 



where the expressions for the elements of J are 



Ja 



pa " 1 



" ^2 



dada 



Jab — Jba 



Jaa — Jaa 
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dhda 
dada 



{a + bxiY 



E 



+ (<'-i)Ei 



xjexp {ti 



i=l 



exp {ti 



^ (1 



x?exp {2ti 



^ (1 - exp {ti) 



'w = Ei 



1=1 



{a + bxiY 

Xi exp{ti) 
exp {ti 



+ 



{a — 1) exp (tj) {a — 1) xf exp {2ti 



exp {ti 



^ (1 - exp (ti) y 
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q2 '"■ 



dpda 

1=1 

_ d'^ ^ xj (g-l) ^ a;fexp(tj) (a - 1) ^ x|exp(2ti) 

dbdb ^ ^ ^ (a + 4 ^ 1 - exp iU) 4 ^ (1 _ exp (t,) ; 
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x?exp (ti) 

Jba = Jab = o o, ^ (^) = o / ^ 77" 

aaao 2 1 — exp [ti 



1=1 

i=l 

Jaa = ^£ {0) = n^' {a + P)- n^' {a) 

Ja0 = Jpa = ^WK-^ {0) = nip' (« + P) 

opoa 

where ■«/''(•) is the trigamma function. 

Under conditions that are fulfilled for parameters in the interior of the parameter space but 
not on the boundary, asymptotically 



V^(d-e)r^N^ (o, , 



where / (0) is the expected information matrix. This asymptotic behavior is valid if / {6) is 
replaced by J{9), i.e., the observed information matrix evaluated at 0. 

For constructing tests of hypothesis and confidence region wc can use from this result. An 
asymptotic confidence interval with confidence level 1 — 7 for each parameter ^j, i = 1, 2, 3, 4, is 
given by 

(§i - z^/2v7^ , + 2^/2 Vj^) , 

where is the ith diagonal element of J{0) and is the upper 7/2 point of standard normal 
distribution. 



5 Application of BLFR to real data set 

In this section, we provide a data analysis to see how the new model works in practice. This 
data set is given by Aarset (1987) and consists of times to first failure of fifty devices. The data 
is given by 

0.1, 0.2, 1, 1, 1, 1, 1, 2, 3, 6, 7, 11, 12, 18, 18, 18, 18, 18, 21, 32, 36, 40, 45, 46, 47, 50, 
55, 60, 63, 63, 67, 67, 67, 67, 72, 75, 79, 82, 82, 83, 84, 84, 84, 85, 85, 85, 85, 85, 86, 86. 
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In this section we fit BLFR, GLFR, LFR, GR, GE, Rayleigh and exponential models to the 
above data set. We use the maximum likelihood method to estimate the model parameters and 
calculate the standard errors of the MLE's, respectively. The MLEs of the parameters (with std.), 
the maximized log-likelihood, the Kolmogorov-Smirnov statistic with its respective ^^-value, the 
AIC (Akaike Information Criterion), AICC and BIG (Bayesian Information Criterion) for the 
BLFR, GLFR, LFR, GR, GE, Rayleigh and exponential models are given in Table [TJ 

We can perform formal goodness-of-fit tests in order to verify which distribution fits better 
to the first data. We apply the Anderson-Darling (AD) and Gramrvon Mises (CM) tests. In 
general, the smaller the values of AD and CM, the better the fit to the data. For this data set, 
the values of AD and CM statistics for fitted distributions are given in in Table [TJ 

The empirical scaled TTT transform (Aarset, 1987) can be used to identify the shape of the 
hazard function. The scaled TTT transform is convex (concave) if the hazard rate is decreas- 
ing (increasing), and for bathtub (unimodal) hazard rates, the scaled TTT transform is first 
convex (concave) and then concave (convex). The TTT plot for this data in Fig. [4] shows a 
bathtub-shaped hazard rate function and, therefore, indicates the appropriateness of the BLFR 
distribution to fit this data. The empirical distribution versus the fitted cumulative distribution 
functions of BLFR, GLFR, LFR, GR, GE, Rayleigh and exponential distributions are displayed 
in Fig. m 

The results for this data set show that the BLFR distribution yields the best fit among 
the GLFR, LFR, GR, GE, Rayleigh and exponential distributions. For this data, the K-S 
test statistic takes the smallest value with the largest value of its respective p-value for BLFR 
distribution. Also this conclusion is confirmed from the values of the AIG, AICC and BIG for 
the fitted models given in Table [U and the plots of the densities and cumulative distribution 
functions in Fig. [H 

Using the likelihood ratio (LR) test, we test the null hypothesis Hq: GLFR versus the 
alternative hypothesis HI: BLFR, or equivalently, Hq: b = versus Hi: b ^ 0. The value of the 
LR test statistic and the corresponding p-value are 3.4 and 0.019, respectively. Therefore, the 
null hypothesis (GLFR model) is rejected in favor of the alternative hypothesis (BLFR model) 
for a significance level > 0.019. For test the null hypothesis Hq: LFR versus the alternative 
hypothesis Hi: BLFR, or equivalently, Hq: (a, b) = (1, 1) versus Hi: (a, b) ^ (1, 1), the value of 
the LR test statistic is 15.3 (p-value = 0.00047), which includes that the null hypothesis (LFR 
model) is rejected in favor of the alternative hypothesis (BLFR model) for any significance 
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level. We also test the null hypothesis Hq: GR versus the alternative hypothesis Hi: BLFR, 
or equivalently, Hq: (a, 6) = (1, 1) versus Hi: (a, 6) / (1, 1). The value of the LR test statistic 
is 8.3 (p- value = 0.0158), which includes that the null hypothesis (GR model) is rejected in 
favor of the alternative hypothesis (BLFR model) for a significance level > 0.0158. For test 
the null hypothesis Hq: GE versus the alternative hypothesis Hi: BLFR, or equivalently, Hq: 
{f3,b) = (1,1) versus Hi: {/3,b) ^ (1,1), the value of the LR test statistic is 19.2 (p-value = 
6e-05), which includes that the null hypothesis (GR model) is rejected in favor of the alternative 
hypothesis (BLFR model) for any significance level. 

6 Simulations 

This section provides the results of simulation study, simulations have been performed in order to 
investigate the proposed estimator of a, /3,a and b of the proposed MLE method. We generated 
10000 samples of size n = 30, 50, 100 and 200 from the BLFR distribution for each one of 
the six set of values of (q, /3, a, b) . We assess the accuracy of the approximation of the standard 
error of the MLEs determined though the Fisher information matrix. The approximate values 
of se{a), se(/3), se{a) and se{b) are computed. The results for the BLFR distribution is shown 
in Table [H which indicate the following results: (i) convergence has been achieved in all cases 
and this emphasizes the numerical stability of the MLE method, (ii) The differences between 
the average estimates and the true values are almost small, (iii) These results suggest that the 
MLE estimates have performed consistently, (iv) The standard errors of the MLEs decrease 
when the sample size increases. 

7 Conclusion 

We define a new model, called the BLFR distributions, which generalizes the LFR and GLFR 
distributions. The BLFR distributions contain the GLFR, LFR,GR, GE, Rayleigh and expo- 
nential distributions as special cases. The BLFR distribution present hazard functions with a 
very flexible behavior. We obtain closed form expressions for the moments. Maximum likeli- 
hood estimation is discussed. Finally, we fitted BLFR distribution to a real data set to show 
the potential of the new proposed class. 
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Table 1: MLEs (STDs) of the fitted distribution, K-S, p-values, -21og(L), AIG, AICC, BIG, 
AD and GM corresponds to times to first failure. 





Distribution 




BLFR 


GLFR 


LFR 


GR 


GE 


Rayleigh 


Exp. 


a 


O.oo47 


A COOT 

0.5327 




0.3520 


0.7798 






(s.e.) 


(0.1432) 


(0.1145) 




(0.0559) 


(0.1351) 






6 


0.1243 














(s.e.) 


(0.0722) 














a 


0.0172 


0.0038 


0.0136 




0.0187 




0.02189 


(s.e.) 


(0.0354) 


(0.0030) 


(0.0038) 




(0.00363) 




(0.00309) 


/3 


0.0348 


0.00031 


0.00024 


0.00031 




0.00064 




(s.e.) 


(0.0025) 


(0.00008) 


(0.0001) 


(0.00008) 




(0.00009) 




-21og L 


460.8 


466.3 


476.1 


469.1 


480.0 


528.1 


482.2 


AIG 


468.8 


472.3 


480.1 


473.1 


484.0 


530.1 


484.2 


AICC 


469.6 


472.8 


480.4 


473.4 


484.2 


530.2 


484.3 


BIG 


476.4 


478.0 


484.0 


477.0 


487.8 


532.0 


486.1 


K-S 


0.1554 


0.1830 


0.1768 


0.2009 


0.2042 


0.2621 


0.1911 


P-value 


0.1786 


0.0703 


0.0877 


0.0353 


0.0309 


0.0021 


0.0519 


AD 


1.749 


2.4890 


4.0346 


3.0923 


3.2530 


13.3205 


3.6505 


GM 


0.3574 


0.4959 


0.5443 


0.6111 


0.6472 


0.8728 


0.6006 
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Figure 1: Plots of pdf of the BLFR distribution for selected parameters. 
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a = 3,(! = 0.1 
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Figure 2: Plots of hazard rate function of tlie BLFR distribution for selected parameters. 
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Table 2: The averages of the 10000 MLE's and mean of the simulated standard errors for BLFR 
distribution, 







AE 


SD 


n 


[a. 0. a, b) 


a 


13 


a 


b 


sd{a) 


sd{l3) 


sd{d) 


sd{b) 


30 


(.5,.5,1,1) 


0.534 


0.537 


2.453 


2.363 


0.449 


0.606 


4.346 


3.454 




(.5,.5,1,2) 


0.528 


0.671 


1.896 


4.485 


0.290 


0.875 


3.191 


5.528 




(.5,.5,3,1) 


.4915 


.7123 


3.321 


3.082 


0.145 


0.571 


5.336 


5.236 




(1,2,1,3) 


1.158 


2.058 


1.977 


5.266 


0.534 


1.415 


3.576 


6.547 




(3,2,1,1) 


4.995 


4.043 


1.275 


2.623 


8.496 


2.761 


1.754 


7.102 




(3,3,3,3) 


3.251 


4.035 


3.249 


3.568 


1.627 


2.054 


3.520 


3.978 


50 


(.5,.5,1,1) 


0.497 


0.714 


1.623 


1.705 


0.151 


0.843 


2.505 


2.153 




(.5,.5,1,2) 


0.504 


0.709 


1.798 


3.720 


0.170 


0.934 


2.233 


4.437 




(.5,.5,3,1) 


.491 


0.730 


3.579 


1.937 


.1165 


.5404 


6.525 


2.658 




(1,2,1,3) 


1.079 


2.053 


1.658 


5.232 


.3852 


1.392 


2.272 


6.470 




(3,2,1,1) 


4.308 


4.237 


1.222 


2.635 


6.246 


2.982 


1.654 


7.820 




(3,3,3,3) 


3.086 


3.874 


3.131 


3.859 


1.307 


1.923 


3.073 


4.716 


100 


(.5,.5,1,1) 


.4916 


.7806 


1.683 


1.293 


.1087 


.9284 


2.698 


1.179 




(.5,.5,1,2) 


.4892 


.7829 


3.597 


3.547 


3.277 


.2104 


5.221 


3.875 




(.5, .5,3,1) 


.4923 


.7513 


3.741 


1.275 


.0815 


.5388 


5.177 


1.556 




(1,2,1,3) 


1.023 


2.169 


1.396 


5.047 


.2841 


1.443 


1.470 


6.017 




(3,2,1,1) 


3.385 


4.122 


1.079 


2.169 


2.669 


2.977 


1.339 


4.954 




(3,3,3,3) 


2.998 


3.779 


3.091 


3.612 


.8535 


1.973 


2.374 


4.357 


200 


(.5, .5, 1,1) 


.4897 


.7999 


1.465 


1.121 


.0739 


.9091 


2.022 


.8251 




(.5, .5, 1,2) 


.4881 


.7773 


1.642 


2.495 


.0845 


.9556 


2.227 


1.902 




(.5,.5,3,1) 


.4936 


.7369 


3.762 


1.036 


.0567 


.4984 


4.427 


1.038 




(1,2,1,3) 


1.005 


2.062 


1.264 


4.522 


.1967 


1.122 


.8104 


4.719 




(3,2,1,1) 


3.070 


4.021 


.9440 


1.987 


1.311 


2.944 


.9673 


3.587 




(3,3,3,3) 


2.999 


3.766 


3.053 


3.626 


.6332 


2.025 


1.940 


3.962 
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X 

Figure 3: The histogram of a generated data set with size 100 and the exact BLFR density (left) 
and the empirical distribution function and exact distribution function (right). 
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Figure 4: Plots of the empirical scaled TTT transform (left), and the empirical distribution 
versus the fitted cumulative distribution functions (right). 
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